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Major  Omeral  Androw  T.  KcNanara 
The  Quarteraaater  General 
Washington  25 t  D,  C, 

Dear  General  NcKanara: 

This  report,  "Predictive  Methods  in  Topographic  Analysis: 

I.  Relief  ,:  Slope,  and  Dissection  on  Inch-to-the-MUe  Maps  in  the 
United  States , "  is  the  first  of  a  ■  series  which  will  supply  the 
Quartemster  Corps  and  the  Anay  with  quantitative  terrain  infor- 
nation  and  with  nethods  for  predicting  the  topographic  cbTracteris- 
tics  of  snail  areas  of  the  earth’s  surface.  This  research  will 
result  in  a  nunerloal  intern  for  describing  the  geometric  dimensions 
of  parts  of  the  earth's  surface.  It  will  also  supply  nsthods  for 
producing  the  infomation  required  by  the  system. 

At  the  present  time,  descriptions  of  terrain  are  expressed  in 
qualitative  terms  which  everyone  interprets  in  the  light  of  his  owi 
experience.  This  report  demonstrates  that  it  is  possible  to  des¬ 
ert  be  terrain  with  nuabers  which  have  the  sane  asaning  for  everyone. 
Quantitative  descriptions  of  terrain  can  supply  information  needed 
in  tJie  writing  of  criteria  for  designing  specific  equipment  (for 
6X&Bipl0f  Vfihl.cl68)i  8nd  In  th®  fcxmolatiion  of  logistical  v)d  tActl* 
cal  plans* 
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foreword 


Pt«15«  rawrlcal  definitions  of  terrsln 
In  order  that  terrain  Infoimtion  aay  be  aatheaatlcaUy 
probleas  as  the  design  of  nilitary  eq«lp«nt  and  logisUcal 
present  terrain  analysis  study  has  been  undertaken  to  attain  this  objec 

tiva. 

Prior  to  establishing  a  nuMrlcal  aystea  for  describing  land^i^s, 
the  fundamental  relationships  existing  among  terrain  T 

understood.  Also,  methods  must  be  developed  for  obtainii^ 
data  quickly,  not  only  ftxm  existing  maps  but  for  areas  ^eie 
mapp^g  is  not  available.  A  series  of  reports  is  proposed  to  furnish 

such  information.  vj 

This  report,  the  first  of  the  series,  describes 
of  six  geomorphic  factors  and  introduces  formulae  for  predicting  throe 

of  them. 
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Abstract 


Based  on  a  200-case  sample  of  United  States  topographic  maps,  a 
study  haiiSUn  made  of  the  inter-relationships  of  six  terrain  factors: 
highest  elevation,  lowest  elevation,  relief,  contour  counts,  slope 
direction  changes,  and  hilltops.  Tto  rank  method  of  correlation  was  used 
for  this  analysis.  From  the  preliminary  correlatiwis ,  three  factors 
(relief,  contour  coimts,  and  slope  direcs  ^on  changes)  were  found  to  be 
most  significantly  related.  Further  worl  with  these  three  factors 
resulted  in  a  8eid.es  of  regression  equat^  ns,  both  siqple  and  multiple, 
to  be  used  for  predicting  any  of  these  t-  elements  idien  pertinent 
data  for  it  are  lacking. 

The  perfoimance  of  two  multiple  regression  equations  to  predict 
average  slope  from  relief  and  from  slope  direction  changes,  were  mapped. 
It  was  seen  that  average  slope  can  bo  predicted  within  an  error  of  ^ 
slope  for  three-fourths  of  the  United  States,  The  best  predictions  were 
made  for  the  Coastal  Plain  and  Interior  lowlands,  areas  of  low  rolief. 

The  poorest  predictions  were  in  the  Appalachians,  parts  of  the  Horthem 
R^kies,  Cascades,  Colors^do  Plateau,  and  Osar^s. 


PREDICTIVE  METHCOS  IN  TOPOGRAPHIC  ANALYSIS:  I*  RBLLEf,  SIOPE, 

and  dissection  on  incH“T0-thb-mii£  maps  in  the  united  states 


1,  Pfflnooae  and  scope 

The  Army  operates  on  the  earth's  surface  and  therefore  conslderstto 
description  of  this  surface  to  be  of  vital  concern.  Heretof^, 
has  been  described  in  qualitative  terms  -  steep, 

which  are  subject  to  misinterpretation  and  are  inadequate  ^ or  ^  g^ses. 
There  is  a  need.  then,  for  explicit  description  of  ^ 

surface  of  the  earth  can  be  measured,  quantitative  descriptions  of  this 
surface  are  possible  and  will  provide  much  needed  information. 

Quantified  terrain  data  can  be  applied  to  the  develoiwnt  ^ 
for  the  design  and  testing  of  military  equipment.  Also,  this  knoiaj^ 
can  provide  a  uniform  method  for  evaluating  field  perfOTnance  of 
ment^ ^tems  with  respect  to  terrain.  Moreover, 

of  su^ factors  as  "relief  and  "average  slope"  ^  is 

studies  of  the  environment  such  as  meteorolo^r 

known,  for  instance,  that  the  roughness  of  the  earth  ® 

the  speed  and  direction  of  wind  and  that  there  is  a  relationship  between 

siae  of  soil  particle  and  angle  of  slope. 

To  fulfill  these  many  needs,  a  numerical  system  for  describing 
scapes  should  be  devised.  This  aystem  should  be  simple,  yet 
enou^  to  carry  all  pertinent  information.  Adaptability  to  “®^® 
cessing  of  these  data,  along  with  those  on 

enviroiuBental  data,  is  also  desirable.  Such  a  aijpq>le  and  meaningful 
system,  or  code,  must  be  based  upon  a  know^e  of  the 
among  Individual  characteristics  of  the  earth' s  surface,  ^e  task  of 
supplying  terrain  information  and  designing  a  ^ 

simple  if  the  earth  were  mapped  entirely  at  an 

Witt  a  20-foot  contour  Interval.  Unfortunately,  the  o^  Mps  available 
for  much  of  the  world  are  at  a  scale  of  1:1,000,000  (15.8 
inch)  with  a  1.000-foot  contour  interval.  Thus,  the  lack  of^tti^ 
topographic  maps  complicates  the  task  of  collecting  basic  terrain  data. 

The  research  here  desc^d-bed  was  undertaken  to  discover  Interrelation¬ 
ships  among  terrain  elements  and  to  demonstrate  how  prediction  equations, 
which  will  help  to  overcome  the  limitations  of  inadequate  map  coverage, 
may  be  constructed,  applied,  and  tested.  In 

development,  the  prediction  methods  have  some  militaiy  uti^y  ®® 

more  S  is  done.  Vm  applicability  of  the  findings  will  become  more 

precise. 


.  ,  ■r.nnvtx''  '  ''  - 


This  report,  the  first  of  a  series  presenting  the  residts  of  a  ooo-> 
tinuing  terrain  analysis  study,  will  deal  only  with  interrelationships  of 
terrain  factors  and  prediction  Methods,  but  the  ultimate  objective  of  a 
code  was  jn  wind,  as  the  research  proceeded.  The  paper  is  further  Ujdted 
to  a  consideration  of  those  topographic  characteristics  which  can  be  col¬ 
lectively  called  surface  geometry. 

2.  PEfy^.9u?,  Uly?g^itttjL9ng 

That  order  exists  on  the  earth's  surface  was  first  demonstrated  by 
John  Playfair,  a  professor  of  matheinatics  at  the  University  of  Edinburgh. 

In  1802,  Playfair  published  a  commentary  on  Hutton's  "Theory  of  the  Barth* 
in  which  he  (Playfair)  stated:  "Every  river  appears  to  consist  of  a  aain 
trunk,  fed  from  a  variety  of  branches,  each  running  in  a  valley  propor¬ 
tioned  to  its  size,  and  all  of  then  toinether  forming  a  system  of  vallles 
( sic }i  communicating  with  one  another,  and  having  such  a  nice  adjustment 
of  their  declivities  that  none  of  ^em  Join  the  principal  valley,  either 
on  too  high  or  too  low  a  level.  *w)  This  observation  is  known  as 
"Playfair's  law"  of  concordant  stream  junctions  and,  although  ejqiressed 
qualitatively,  it  has  been  the  basis  for  recent  hydrological  investlgaticne. 

A  quantitative  interpretation  of  Playfair's  I«v  was  iindertaken  by  the 
hydrologist,  R,  E.  Horton,  In  19*+5i  Horton  showed  that  stream  channels 
develop  according  to  definite  hydrophysical  laws.  Thius  he  was  able  to 
express  the  conposition  of  the  stream  system  of  a  drainage  basin  quanti¬ 
tatively  in  terms  of  stream  order,  drainage  density,  bifurcation  ratio, 
and  streaav-length  ratio.  0)  Since  that  time,  other  hydrologists  have 
elaborated  vpon  and  refined  the  basic  conc^ts  which  Horton  descriLbed, 

These  first  quantitative  studies  were  concerned  with  stream  systems. 
However,  if  stream  systems  have  oz^er,  their  interfluves,  ^ich  have  been 
shaped  by  them,  must  also  exliibit  order.  Recently  geologists  and  geog¬ 
raphers  have  turned  to  a  quantitative  interpretation  of  other  elements  of 
the  earth's  surface.  In  1950,  Strahler  used  statistical  methods  to  test 
alternative  hypotheses  in  geology.  (?)  Other  studies  by  Strahler  and  his 
associates*  have  made  it  abundantly  clear  that  the  apparent  disorder, 
which  had  previously  been  seen  in  the  geometry  of  the  earth's  surface,  is, 
in  reality,  order. 

From  the  above  statements  it  should  follow  that  randomly  chosen  units 
of  the  earth's  surface  should  have  mathematical  relationships  among  their 
separate  dimensions ^  Also,  it  is  likely  that  information  about  large  areas 
gathered  from  small-scale  maps ,  would  have  a  relationship  to  the  smaller  , 
areas  contained  within  them.  If  this  is  true,  the  possibility  exists  that 
by  using  regression  relationships,  'predictions  can  be  made  about  small 
areas  which  have  not  been  mapped  at  topographic  scale. 


*  Working  under  Office  of  Nav?!  Research  Contract  N6  ONR  271. 
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So  far  as  is  knowi,  only  three  probing  attempts  have  been  made  to 
discover  whether  it  is  possible  to  develop  regression  equations  for  pre¬ 
dicting  terrain  characteristics  on  randomly  chosen  parts  of  the  oaHsh's 
surface.  Peltier,  in  195^,  found  that  average  slc^  is  oapirtcally 
related  to,  and  varies  directly  with,  the  average  local  relief  frosi  one- 
square-nile  areas* (5)  The  following  year,  another  investigation  showed 
that  regression  eqxiations  ior  predicting  the  relief  on  a  squaire  nile 
became  Incx^singly  more  accurate  as  a  set  of  terrain  sanples  was  divided 
and  subdivided  on  the  basis  of  natural  vegetation  type  and  geologic  bed- 
rock.'iO'  A  third  study,  one  towards  developing  prediction  equations,  was 
made  in  1956  when  it  was  demonstrated  that  average  slope  on  rather  small 
areas  (9  square  kilcsneters)  could  be  predicted  through  the  use  of  data 
gathered  from  maps  with  scales  of  1:250,000  and  1:1,000,000.(9/  These 
three  studies  were  small  in  terms  of  the  amoxint  of  data  involved  and 
their  only  purpose  was  to  investigate  the  possibility  of  developing  sta¬ 
tistical  methods  for  predicting  relief  or  average  slope  of  randomly 
chosen  small  areas. 

The  results  of  the  works  cited  above  were  promising,  and  a  much  more 
ambitious  project,  involving  more  data  and  much  more  extensive  analysis 
was  undertaken  by  the  present  authors.  Means,  modes,  medians,  and 
standard  deviations  of  the,  te^crain  factors  to  be  treated  in  this  latest 
study  have  already  been  published  .(11/ 

3.  Present  otii^ntitative  st\>dy 

a.  Basic  data 

The  sample  selected  for  study  consi.sts  of  a  group  of  204  US  topo¬ 
graphic  maps  chosen  at  random  from  the  Amy  Map  Service  index;  these  maps 
are  at  a  scale  of  1:62,500  or  1^50 *000 ,  These  sheets  represent  every 
major  physiographic  province  of  the  United  States  and  include  a  wide 
variety  of  terrain  features  (Pig.  l). 

From  this  sanple,  data  were  gathered  for  six  terrain  elements.  These 
were : 

By  area: 


!•  Highest  elevation  within  a  sample  unit. 

2,  Lowest  elevation  ^thin  a  sanple  unit. 

3,  (difference  between  the  highest  and  lowest  elevations 
of  the  sample  \init)  within  a  sanple  unit. 

4,  Number  of  hilltops  (individual  heights  of  land  as  represented 
by  a  closed  coiitoxir)  within  a  sample  unit. 


3 


Figure  1.  Distribution  of  the  204  topographic  saps  of  the  United  States  which  provide  a  reprasen- 
tative  terrain  saople. 


fiy  traverse; 

5»  Kuinber  of  20-foot  contour  crossings  (contcx^j*  counts)  along 
a  inverse  within  the  saaple  unit,  used  as  a  neasuz*e  of 
average  slope. 

6.  Number  of  valleys  and  divides  encountered  on  the  sai.«  tra¬ 
verse,  used  as  a  measxire  of  dissection  and  collectively 
referred  to  as  slope  direction  changes. 

The  smallest  sample  unit  is  a  circle  with  an  area  of  ^/l6  of  a  square 
mile  I  the  mid-point  of  hhich  is  the  center  of  the  topographic  sheet. 
Successive  doublings  of  this  unit  produced  sample  circles  with  areas  of 
5/8 1  li.  2^,  5,  10,  20,  40,  80,  am  l60*  square  miles,  resulting  in  10 
sau^le  sizes.  Two  diameters  were  drawn  in  the  cardinal  directions 
through  the  mid-point  of  the  concentric  circles.  These  lines  were  used 

;  "  as  traverses  along  >rtiich  contours  and  slop>e  direction  changes  were  counted. 

I  The  combined  lengths  of  the  two  diameters  for  each  circle  are  1.2,  1.8, 

I  2.5,  3.6,  5.1,  7.1,  10.1,  14,3,  20.2  and  28.5  miles  (Fig.  2). 

Thus  there  ware  ten  area  sizes  per  sheet,  with  a  total  of  2,040  indi¬ 
vidual  sanqjle  xinits  to  be  considered.  For  each  sample  unit  there  were 
I  gathered  the  previously  enumerated  information  on  the  six  terrain  ele- 

I  ments,  making  a  total  of  12,240  pieces  of  data  for  analysis, 

I  b.  Preliminary  analysis 

first  analysis  of  these  data  revealed  an  orderly  pi'ogression 
a  in  the  change  of  value  of  individual  measures'  with  change  in  siie  of  area, 

I  The  data  also  exhibited  a  strong  skewness  which  all  attempts  failed  to 

I  normalize.  However,  if  several  statistical  distributions  refuse  to  be 

I  nonnalized ,  they  can  be  transformed  into  identical  rectangular  distribu- 

5  tlons  prior  to  analysis.  This  will  assure  that  PJ.t.reme  values  within  the 

^  population  do  not  exert  an  undue  influence  on  ai.'V  regression  equation 

which  may  be  developed, 

‘t 

j  A  quick  analysis  of  twenty  cases  drawn  from  the  total  samples  indl- 

I  cated  the  feasibility  of  using  rank  correlation  methods  for  ease  of 

interpreting  data.  The  relief  values  of  these  cases  were  ranked  for 
{  each  size  of  area,  i.e, ,  arranged  in  ascending  order,  of  magnitude  and 

I  assigned  a  value  from  1  to  20,  v/ith  the  lowest  value  of  relief  becoming 

i  1  and  the  highest  beccaning  20.  These  rank  values  for  each  size  of  area 

‘  were  then  cori'elated  with  those  of  every  other  size  area.  The  coefficients 

of  rank  correlation  showed  an  orderly  change  of  value  with  change  in  area 
i  of  samples,  thus  indicating  the  degree  of  success  to  be  expected  if  relief 

I  of  a  large  area  is  used  to  predict  relief  of  a  small  area  which  lies  at 

its  center. 


*  A  few  of  the  l60-square-mile  samples  had  to  be  arranged  as  elipses  to 
stay  within  the  limits  of  the  topographic  map. 
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MESSENA.  NEW  YORK 


CONTOUR  C0UNTS(4i 

'slope  direction  changes 


TRAVERSE 


5/8  sq.  ni. 

2.5  aq«  Bi. 
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Hi^iest  Elevation  •  feet 
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Lowest  Elevation  >  feet 
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12 

Figure  2.  Portion  of  O.S.G.S,  topographic  nap  illustrating  data- 
gathering  technique.  Highest  and  lowest  elevati<xis  reed  to  the 
nearest  contour.  When  using  a  topographic  nap  for  the  gathering  of 
data,  interpretations  are  node  for  elevations  both  above  highest 
contour  and  belov  lowest  contour. 
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It  was  also  shovm  that  th«  values  of  the  coefficients  of  correlation 
increased  in  an  orderly  fashion  as  the  sire  of  the  independent  variables 
approached  the  size  of  the  dependent  variable.  This  orderliness  was  so 
pronounced  that  for  the  purpose  of  futiw  study,  which  would  be  done  by 
rank  correlation  nethods,  it  %ras  deemed  safe  to  eliminate  half  of  the 
saaple  according  to  size  of  areas  and  concentrate  only  on  the  5/8,  2^, 
10,  40,  and  l60-square-mjle  samples. 

Relaticffiship  among  terrain  elements 

a.  Correlation  bv  the  rank  method 


Rank  correlation  does  not  offer  the  only  possible  method  for 
analyzing  the  data.  Graphical  curvilinear  methods,  for  instance,  would 
be  applicable,  and  in  some  cases  give  more  accurate  results.  However, 
such  an  approach  is  labozdous  and,  when  more  than  one  independent  variable 
is  used,  extremely  slow.  By  adc^ting  the  rank  corx*elation  method,  it  is 
possible  to  put  the  data  on  punch  cards  and  do  much  of  the  confutation 
on  machines,  thereby  gi^atly  shortening  the  task.  Strictly  speaking,  for 

this  work,  the  Pearsonian  equation  for  correlation  rr  =  -mJimSLmmmm  was 
used  with  rank  values  substituted  for  absolute  values  Instead  of  the  usual 


equation  for  computing  rank  correlations  rj*  =  1  - 


6  Ld2 

N(n2  -  i) 


The 


Pearsonian  formula  is  more  readily  applied  to  the  confutations  furnished 
by  the  EAM  machines.  When  200  paired  values  are  used  for  a  correlation, 
tee  results  are  so  nearly  the  same  that  accuracy  is  not  sacrificed, 

(Rank  correlation  of  contour  counts  and  slope  direction  changes  on  5/6 
square  miles:  Pearsonian  formula;  r  =  .5195;  rank  formula:  rr  “  .5228,) 


Simple  correlations  reveal  the  degree  of  interrelationship  among  indi¬ 
vidual  factors  and  indicate  which  ones  can  be  used  to  predict  others.  The 
degree  of  interrelationship  is  indicated  by  the  end  product  of  correlations, 
the  value  rr.  This  measxaro  notes  the  significance  between  elements  and 
can  be  evaluated  from  the  charts  or  tables  of  reliability  found  in  any 
standard  statistics  book. 


For  this  studj'  of  200  cases  (for  ease  of  statistical  manipulation, 
four,  of  the  original  204  samples,  wore  discarded)  an  rr  value  of  .1235  is 
required  for  the  5^  level  of  .significance* ,  according  to  the  F  table  of 
significanctsi, However,  for  one  factor  to  be  a  good  predictor  of 
another,  the  value. must  be  considerably  higher,  A  useful  indicator 
of  prediction  accuracy  is  the  rr^  which  will  tell  what  percentage  of  total 
variance  between  two  elements  has  been  accounted  for  by  the  coiTelation, 
Thus  an  rr  value  of  ,1235»  although  significant  in  that  it  probably  did 
not  occur  by  chance,  accounts  for  only  2^  of  the  total  variance 


*  This  value  of  rr  would  occur  only  once  in  20  times  by  chance 


TABU  X 


Slope 

Ureetioa 

Chaagaa 

Highaat 

ElevaUon 

Lowaat 

Hlavation 

Relief 

HiUtopa 

Coatour 

Counta 

S/8  Eeuara  Mlea 

Bigh  Elevetion 

Low  Hl'svB’-lon 
Relief 

Hilltopa 

Contour  Couute 
Slope  Dir  Cheocea 

■  9799 

5507 

.h2h7 

mmmm 

.0526* 

.0366* 

.2213 

•  V  W4a 

.505*^ 

.3663 

.9‘*72 

.2935 

.0782* 

.OA39* 

.3957 

.*•078 

.5193 

2  1/2  Sauara  lUlaa 

Ugh  Elevation 
low  Elevation 
Elilief 

Hilltopa 

Contour  Counta 
Slope  Dir  Cbaafaa 

.9619 

.6277 

.h5^ 

.0716* 

.0538* 

.2395 

5235 

.^60*^ 

•9291 

.3559 

.1106* 

.0677* 

.3752 

.52^ 

.5567 

10  Sauara 

Nilea 

High  Elevation 

X^  Elevation 
Relief 

HiUtopa 

Contour  Counta 
Slope  Ur  Changea 

M  » 

.9‘«02 

mmmm 

.7009 

.k755 

.073^* 

.0308# 

.1967 

•  e»aaa» 

/ 

.5357 

.3176 

.09U 

.3715 

.0305* 

-.0299* 

.2684 

.5327 

•  5332 

'  mmmm 

40  Square 

Mi  lea 

High  Elevation 

Low  Elevation 
Relief 

Hilltopa 

Contour  Counta 
Slope  Dir  Changea 

.9006 

.7707 

.^969 

.1036* 

.o656« 

.2025 

.5712 

,3lW 

.8632 

.3*»71 

-.0238* 

-.0466* 

.1549 

.5006 

.4651 

l60  Sauara  Kllea 

High  Elevation 

Low  Elevation 
Relief 

Hilltopa 

Contour  Counta 

mmm"^ 

.8619 

.ak6i^ 

•  5215 

.1381# 

.1071* 

.1935 

.6079 

.31^3 

.8267 

.3029 

-.0606* 

-.0555* 

.0671* 

.4239 

.4153 

Slope  Dir  Cnuuices 


*  not  Blgidflcent  tt  at  ^  level  (rr  •  .1235  required  for  •Ignlflcance) 

#  not  elgniflcant  rr  at  1$  level  (rr  -  .1525  required  for  elgnlfioence) 
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COmTOM  COUMItCOMTOU*  COVMTI 


COCrVCftNTS  Of  tANl  COIKIATIO** 


itJH 

ta 


1.12352  ■  .0152)  Vtta 

ui  rr  of  .92  account#  _ti  t«min 

purpososi 

Wh.n  ««eh  Urriln  ^«°^rt5Uc  ii 

*  ^  ^  t«KLa  I  as  wall  as  other  considerations ,  permitted 

k  close  study  of  Table  i,  as  w^  m  *l«vatioo,  and  hilltops  trem 
the  ellMlnatlon  of  hig^st  louest  elevation  have  high  correla- 

further  study.  relief.  (Klevations  correlate  sell 

tlons  ofOy  with  is  ihe  difference  between 

with  relief,  )  Althou^  •hilltops"  (defined 

highest  elevation  and  l^^jt  tlMS^ont<«r)  la  an  i>,>ortant 

as  any  siuMlt  repwaanted  ^  sliiriJiated  froa  further  study, 

elenent  of  the  landsoye,  this  -Mialtlve  to  differences  in  con- 

s:;:  -  :rT,^tro?  ^c. 

Ute  wall  with  each  other  STsample  alaea  being  cen- 

relatlonshlps  fj?  ^  fTeble  H)  From  this  table  It  is  noted 

sidered  are  widely  different  (Table  W.  w  relatively 

that  th.  rj.  TOlu..  ff.‘!2trSn.M^T.ln  riMBUt  ot  MOttMr  .U.  of 
little  when  correlated  with  of  relief  and  contour  counts 

area.  For  Instance,  an  rr  value  of  ^72  of  on  5/8 

The  change  In  the  values  of  construction 

change  in  size  of  area  is  «o  ^'^"bably  have  been,^ 

of  graphs  which  show  what  u-od.  (Coefficients  of 

slzM  of  area  other  than  t^  five  ^%tied  for  relief, 

rank  correlation,  in  relation  to  area,  have  l^en  plo^;^^ 

contour  counts,  and  sl^  iell^wlief^  10  square  niles 

exai^le,  one  adght  ^nt  Figure  3  selwit  the 

.  correlates  Tilth  relief  on  80  .  r«Tii-«»sants  80  square  ailes. 

line  on  the  botton  ^^represents  10  square  ailes. 

Follow  this  line  upwa^  S^vThedetei^ed.  The  correlation 

rSirr«  ?0  to  ..li.f  =n  80  ,<r-«  -^.s 
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1 


is  shout  .94.  Ihs  othw  gruphs  (Figs.  4.  5.  •«»  6)  will  supply 
polstad  Infoimtlon  rsgsrdlng  contour  ecunts  snd  slops  dlrsctlon  chsngss 
in  the  saaa  Mnnar. 

The  correlation  coefficients  for  the  three  factors  on  all 
area  supply  the  basic  information  needed  prior  to  the  making  of  pi^lction 
formilae.  The  tables  of  rr  (Tablee  I  and  II)  indicate  the  sign^cance 
of  one  factor  to  another.  The  numerical  values  of  Tj-  can  be  used  for 
computing  simple  regression  equations  based  on  rank  correlations. 


b.  Predlctionfl 

(1)  Staple  regression  eouatlons.  ^Kression  equations  were 
computed  to  predict  rank  values  for  the  following  combinations  of  terrain 
factors  on  all  sizes  of  area: 


Table 

1. 

Relief  to  relief 

III 

2. 

Contour  county  to  cwitour 
counts 

IV 

3. 

Slope  direction  changes  to 
slope  direction  changes 

V 

4. 

Relief  to  contour  counts 

VI 

5. 

Slope  direction  changes  to 
contour  counts 

VII 

Relief  to  slope  direction  changes  has  br^n  omitted  because 
correlation  values  indicated  that  not  much  of  the  variance  in  one  factor 
can  be  attributed  to  differences  in  the  other. 

The  standard  procedure  of  the  least  squares  methcd  of  nt  was  «ad 
to  estimte  the  regression.  T  =  a  +  b  X.  -T"  /he  prrflcW^  "a  a 
constant,  "b"  the  coefficient  of  regression,  and  X  the  ^e^^ent 
variable.  With  rank  correlations,  the  oonputatlons 

not  complex,  since  "b"  is  identical  to  the  rr  values,  and  a  Is  computed 

from  the  formula:  a  •  (»  *  DU  -  where  -K«  le  the  nuiaber  of  cases, 

and  "rr"  is  the  coefficient  of  rank  correlation. 

The  sanple  size  for  this  study  was  200,  hence  the  regression  formula 

for  any  factor  in  this  study  becomes J  T  »  ^20l)(l  -  r)  +  pj,  x.  Should 

one  wish  to  devise  regression  equations  froa  information  taken  f^  th3 
graphs  of  rr  values  (Figs.  3,  4,  5.  and  6),  this  fomula  would  be  u^^ea. 
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An  «»pl.  of  th.  .ppUctlon  of  “““(J^'^VSiTof  .”lo!iS^o- 

f/8*.;Sn^  ilSTSi^ibft  .i-.  Th.  .^Uon  t  . 

■U«»,  1.691  f.«t  Is  con».rUd  to  s  rsnk  valu.,  17^ 

convort  aetu.1  toIum  of  r.U.f  to  >»nk  valuss  of  ^ 

Ho«^ph,  for  ^ 

for  cont»oui*  countfl  wid  sxopo  dlrdcvicwi  crwis  ®  -• 

9?f  1^5tUul?S^r«  md.  in  the  eq^Uon  and  the  equation  solved. 

T  ■  9.6078  +  C.9044  X  173) 

T  c  9.6078  +  156.46 
T  ■  166.06 

THisultina  a  prediction  of  relief  for  the  5 /8- square-#^  ar«  as  ^ 

S^roS^f'S  m,  si  th.  pr^licuon  is  70  f..t 

more  than  what  actually  occurs* 

Sioflo  regrssslon  equations  to 

ohanees,  and  contour  counts  have  b<wn  dwelopwl  for  the  area  sls«. 
studW  ai«i  it  has  bwn  shovm  that  a^ltional  %TL  ^ 

could  be  dmelopad.  With  these  equations  one  should  ^  abl^  ^ 
fairly  accurate  estiaates  of  these  ela*nts  on  w  »l-“  of  area  uh 
Infonaation  for  only  one  of  the  factors  is  available. 


In  the  preceding  section,  relief,  conW  “"L 

slope  direction  changes  have  been  buUt  indivMually  into  ^ 

through  the  developoant  of  multiple  regression  equations. 

Hultiple  regr.5ssion  is  sisdlar  to  slsple  ”>8res3ion.  but  ^es  UM 
of  JolS  than  onl  factor  for  predicting  another.  J^eojetica^.  « 
y>m  vioastble  to  determine  all  the  iE5>ortant  elenents  that  constitute  an 
individual  piece  of  the  land  surface.  All  these  *^*!*"H 
oluld  then  be  built  Into  .1  matiple  ^slS^ine. 

theoretically  at  Isast.  give  nsar  llntrtbSinTfaltors 

However,  at  thia  stags  of  tsrrain  res^h.  all  ^s  contributing  raooors 

8,re  not  known;  neither  Is  adequate  information  available. 


Figure  7,  For  use  in  the  prediction  equations,  actual  values  of  relief 
on  five  sizes  of  area  are  converted  to  rank  values  of  relief,  or  the 
reverse,  by  means  of  the  nomograph. 


flinire  8  For  use  in  the  prediction  equations,  actual  nuaiber  of  cOTtour 
^  nZ  of  arL  are  concerted  to  rank  number  of  contour 

:ount8,  or  the  reverse,  by  means  of  the  ncanograph. 


•OMMUAMt  TO  COOMOT  MtMMCO  O#  MOOI  DtOgCTlOO  CMMMII  TO  OAHKI 


Figure  9«  ^or  use  in  the  prediction  equations,  actual  number  of  slope 
direction  changes  on  five  sizes  of  area  are  converted  to  rank  number  of 
slope  direction  changes,  or  thie  reverse,  by  means  of  the  nOTOgraph, 


\mm 


=rs;.n4s^  2^  iji^ssuis'/mS  »-. 

Tiirt  t.4ftt  easfls  Haro  p6rfom6d  to  lUttdtrate  this  tochnltju* 
SSchTIdicatS  thoM  alwwnto  «uU  bwt  i^iw.  th«  oi  Klnal  pr«- 
dictions. 


1,  Slope  direction  changes 

Froa  2*  Relief 

3,  Contour  covrots 


Correlation^ 
rr  rr^ 

,3079  *  0^ 

160  sq  miles  .785^  .6172 

.8848  .7829 


The  three  factors  used  in  all  possible  combinations  produce  the  corre¬ 
lations  listed  below. 


Correlation, 


From 


Relief  and  slope 
direction  changes 
Contour  counts  and  sl()pe 
direction  changes 
Contour  counts  and  relief 
Contour  coxints,  slope  di¬ 
rection  char.goB,  and  jrolief. 


160  sq  miles 


.8262  .6827 

.8872  .7871 


.8897 

.8899 


.7916 

.7919 


The  ri.’  s  of  the  simple  regression  show  that  sl<^  direction  changes 

on  Ito  iu...  a,^  prJIctor  of  »Uef  . 

very  coot:  r.ll.f  »s  a  pradletor  of  relief  1»  good.  «id%ontour  coonM 

l.^th.  b«»t  prwllctor  of  rri-lef.  Th.  “rrolation  ocfflcl^ts  rfto. 

rmltlpl.  regi«alon  equations  indicate  ^ 

ll^glor  S’su^?  wJSlq»t[S  oo^taln- 

no  significant  ixaprmmiBnt  over  the  single  r'egreaslon  eq^  i+Mir  \n 
con?^“s  S^e.  Although  "slope 

not  a  good  predictor  of  relief,  when  added  to  either  of  the  other  two 


■o-  >1?; 


factors  It  helps  to  inprova  the  estiaato.  (Relief  by  Ueelf  accounti 
for  only  6l.?2^  variance;  In  co«binatlon  with  slope  direction  chancea, 
68.2^  of  tho  variance  is  accounted  for.) 

ilie  other  series  of  equations  predicta  contour  counts  on 
■ile  sas?)lea  fro«  l60-8quare-eile  inferwition.  The  aystee  used  above  U 

followed : 

Tfi  nr^ict  contour  counts  on  5/6  square  wiles  (sllPlff 

Correlation 

£t _ 


Fro* 


Slope  direction  changes 

2.  Relief 

3.  Contour  counts 


s' 


.4141  .1718 

160  sq  KUes  .72^5  .52^ 

.8955  .8219 


To  predict  contour  cotints  on  5/8  fouare  (unltiole  regrOfgigJll. 


From 


’4,  Relief  and  slope  direction 

changes 

5.  Relief  and  contour  coxints 

6.  Slope  direotion  changes 
and  cpntqor  counts 

7.  Relief,  slope  direction 
changes,  and  contour  county 


1^0  sq  Biles 


Correlation 

a _ 

.8119 

.6591 

.8960 

.8029 

.8967 

.8041 

.8967 

.8041 

estimates 

made  pos- 

sible  by  use  of  multipl.e  regression  aquations.  They  also  show 
count  of  contours  along  a  random  traverse  (Equation  V  c*n,  ^  itsexT, 
yield  nearly  as  good  sji  estiaate  when  used  alone  as  when  used  in  combina¬ 
tion  with  other  data. 


Although  "contour  counts",  which  may  be  converted  directly 
average  slope,  is  one  of  the  most  discriminating  of  terms  deacrtbing 
surface  geometry,  it  is;  at  the  same  time,  the  most  difficult  informa¬ 
tion  to  obtain  fo\'  many  parts  of  the  world.  Multiple  regression  equa¬ 
tions  for  predicting  contour  counts  from  relief  and  slope  direction 
changes  (Equation  4)  on  unit  areas  would  have  practical  use  on  many 
foreign  maps  which  show  ridge  lines,  valleys,  and  spot  elevations,  but 
no  contours.  The  principle  would  be  applicable  also  to  tho  ana^sia  of 
aerial  photographs  from  which  average  slope  readings  aro  more  difficult 
to  obtein  than  relief  values  and  slope  direction  changes. 
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tion  changes  are  presented  in  Table  vill.  atw  enoatlona, 

Vin  are  coapared  with  aUllar  caaea  rf  th«  »iap^  • 

it  la  obvious  how  laoch  the  eatSna tea  have 

wh«n  predicting  contour  counts  on  l60  square  a  3 

SlJectlon  cha^ea.  r,  1.  .9038i  predicting  froc  relief  alone,  rr  la  only 

e72^5  (Table  11  )e 


<yqustaona 


To  tost  t^*•  offectHrenoso  of  the 
theoretical  value  a,  indicated  the  'T* 

unit  areas  were  predicted  by  vising  tho  multiple  regression  equa^wa 
presented  In  Table  VHIe  The  predicted 

nressed  as  an  actual  number  of  contour  counts  wnich  were  then  converveu 
to  a  slope  tangent  by  means  of  the  Wentworth  equation: 

Tan.  =  7?n^^y  ^ 

Tva.  a<iff«rence  between  the  predicted  and  actual  slopes  was  napped  for 
S:  fwo  SIT.  (5/8  and  l60  «,uar.  ail..)  In  t.™»  of 

percent  slope  (Fi.gs.  10  and  11). 

In  ae'ieral  the  equatioiis  perform  well  in  areas  of  low  relief  swh 

u'n”1s"S^for*?2el6S-Jiua«^^S^ 5^-«.na»-alla 
samples,  the  same  areal  pattern  of  errors  occurs. 

However,  if  one  were  to  decide  that  a  prediction  error  of  ^ 

tolerable,  it  could  be  concluded  that  the  equations  us^  to 
^  +  p^rnttaAq  10  and  11  are  adequate  for  predicting  average  slope 

of  200  topographic  maps  adwjuately  jvgiresmts  th.  country  as  a  wh  . 


'■■  Th.  qu.5tion  might  wsU  have  arisen  in  the 

s  rs^u&a!  w 

'Z1  SoSefStfornf  U^lefa^Xtof 

vavanriv-^ivAfi  If  thls  <=ampi  A  doos  Hot  adequately  represent  the  wuriac^ 
«™rtS’of  thrS^i^S^states.  then  the  findings  apply  to  the  eaa,.!. 
g™lSa?ly.  and  only  generally  to  the  rest  of  the  country. 
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Figure  10.  Mapped  differences  between  the  actual  nuiber  of  contour 
counts  on  5/6  square  nile  and  the  predicted  number  (converted  to 
al^) ,  as  derived  from  a  multiple  regression  equation  (!^ble  VIU) , 
indicate  the  areal  di  stributlon  of  the  equation '  s  e  f  f  ec  tiveness . 


HcroMiWM*  or  MuiTifu  nuticMMa 

CQ^iOM  cuuNfi  :  lai;*}  »  aois  iiLai  tiort  b«icTaW|';H#wcii, 


o 

Mapped  differences  between  the  actual  number  of  contour 
counts  on  160  square  miles  and  the  predicted  number  (converted  to 
slope),  as  derived  from  a  multiple  regression  equation  (Table  VIII), 
indicate  the  areal  distribution  of  the  equaticvi'a  effectlvmess. 


-r.  ■' 
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It  la  not  plannod  to  give  a  final  tncl  daflnltlv#  anawar  to  thla  quoiy 
sinco  ths  pofsalblllty  axlsta  that  tha  aan>l*  aay  b«  adaqiiata  In  mm 
raapMta  and  Uiadaquata  in  othars  (l.a..  adaqiiat# 

alopa  on  160  aquara  ailaa,  but  not  adaquata  »dth  ragard  to  raliaf  on  5/0 
sqira  aUaa).  Hoaavar.  it  la  poaaibla  to  obUin  a  flrat  Uprtaaicn  of 

tha  aaapla'a  adaquacy. 

It  My  ba  safaly  aaauaad  that  a  saigila  satisfactorily  rapraaanta  tha 
population  froB  ahich  it  ia  drawn  if  it  aeata  thaaa  taata:  UJ  aaaan- 
tially  tha  aane  Maans  and  atsndard  deviations  ara  ob^^ad  whan  tha  aaapia 
has  baan  randoady  dlvldad  into  two  equal  groups  and  (2)  a  ragrassion 
equation  basad  on  one  half  of  tha  aai^pla  shows  the  saaa  perforranca 
charactariEtics  in  both  halvas. 

Whan  these  tests  "sra  perfor«ed,  it  was  found  that  the  uMns  and 
standard  davtations  o*  tha  two  sub-greops  of  tha  and  the  perxor- 

■ance  of  the  prediction  aquations  ware  slailar.  (Sea  Appendix,;  Th^, 
it  nay  be  said  that  tha  8an}la  is  adequate  in  that  it  represents  a  sta- 
tlsvieal  distribution  of  tha  (population  of  all  of  the  possible  16C-8quare- 
Bila  units  in  the  United  States.  In  other  words,  the  noMgraph 
ihoirtn*  th.  ranking  of  th,  two  hundrel  l60.Muar^ail, 

to  contour  counts  (indicating  average  slope)  would  ohow  little  dlfferar.ee 
froa  a  chart  basad  on  every  l60-aquare-»lle  area  in  the  country. 

The  adequacy  of  this  200-case  sasgile  to  represent  the  tenpin  popula¬ 
tion  of  the  United  States  should  not  be  confused  ulth  the  accuracy  of 
predictions  based  on  it.  It  has  already  been  noted  (Figs.  10  and  11 ; 
that  predictions  of  average  slope  are  not  reliable  in  areas  of  rugn 
relief.  Better  estinates  of  avorigfe  in  these  areas  oiay  require  a 

closer  network  of  samples  and  additional  data  concerning  such  factors  as 
natural  vegetation  and  bedrock. 

5.  Conclusions 

This  repiDrt  has  been  written  at  a  stag©  in  terrain  research  where 
it  is  convenient  to  take  stock  of  the  work  acconqplished  to  date,  noals 
for  future  study,  and  methods  for  achieving  them,  can  be  formulated  on 
the  basis  of  the  findings  published  in  this  reijorti 

Again,  the  orderliness  of  th©  earth’s  surface  has  been  demonstrated, 
this  time  by  the  conversion  of  terrain  data  to  a  rectangular  disl/ribu- 
tion  prior  to  analysis  and  the  presentation  of  the  results  on  gniphs  and 
charts.  More  details  of  the,  pattern  and  order  which  exl.st  among:  land- 
form  features  will  be  revealed  through  extended  analysts  of  thase  sm 
data.  Furthermore,  the  technique  of  conversion  eii?>lqyed  here  prwldes  a 
relatively  quick  and  easy  means  to  manipulate  quantified  terrain  (informa¬ 
tion.  .  '  I  ;  ..  i|  '  , 


f 


With  tho  con^lition  und  5Ub9#^uoDt  of  six  torniln  factors , 

•oiie  of  th<i  relationships  sxlstijijg;  aac^  geoaorphic  •Isasnts  )av«  base 
discoveral.  Moreover,  as  undersUndtiig  of  these  reUUonrtiips  has  eoM 
about.  It  has  been  possible  to  devise  regression  equations  for  pradletlng 
various  terrain  elaaents  where  pertinent  data  aay  be  laeklng.  A  series 
of  regression  equations  has  been  foraulated  for  predicting  relief  and 
average  slope.  Haps  of  the  perforaance  of  these  equations  have  pointed 
out  those  areas  not  likely  to  require  further  stxriy  (insofar  as  relief 

and  average  slope  are  eonoemed)  as  well  as  other  sreas  which  will  require 
further  study. 


The  prediction  equations  presented  here  apply  only  to  the  contlnenUl 
Hnlted  States.  Other  areas  of  the  world  are  also  orderly,  but  their 
order  differs  fro«  that  of  the  Dnited  States  to  the  sane  degrw  as  do 
the  structural  ■ateri.als  and  the  physical  processes  which  ,  have  shaped 
then.  However,  the  sane  techniques  of  sanpllng  and  data  analysis  used 
to  study  the  United  States  can  be  applied  to  other  regions. 
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PRoor  or  swu  ads^ct 


Th«  wi«  divldsd  Into  tvo  froupf  of  100  e«M«  ««eh  bgr  ftXloidng 

di9lnt«r«st«l  co-worinirf  to  ohufflo  th*  daU  card*  and  ditrld*  than  into 
two  jsllot,  labolcd  Oroops  I  and  IX.  tha  aouia  and  auindard  davlationa  of 
iheao  two  groups,  along  with  <1m>s«  of  tha  tot4l  aajf>la,  art  as  follows: 

fotal 

Qrgyp  X  Sism-U 


[Koan  nunbor*  5^2,16  531.3^  536.76 

.. 

Standard  davlatlon  559.67  525*10  536*07 

[Moan  .1132  .1109  .1121 

' 

Standard  dofwlatlon  .1169  .1096  .1121 

*)fean  nuaber  of  20-foot  contour  countar  along  2  dianatara  of  a 
l60-aquAro-ailc  circla. 

The  data  of  Orotjp  I  transfomed  into  a  rectangular  distribution  by 
the  ranking  aethpd  already  demons tratad,  prc^uced  the  following  rogros- 
sion  equation  for  predicting  conto’ir  counts  m  i60  square  miles  trooi 
slope  direction  chan;?.e8  and  relief  on  16C  square  miles:  T  *  -9*0955 
.8137  Xi  +  ,3664  ^2  *-9. 0955^18  a  constant,  "Xi"  is  rank  value  of 

relief  ,  j and  "X2"  is  rank  value  of  slope  direction  changes.  (It  will  be 
seen  t^t  this  aquation  is  nearly  identical  to  the  equation  developed 
for  the  saws'  purpose  from  the  200-case  saiqjle:  T  «  -16.6761  +  .8045  1^1  “*■ 
,3614  X2  (Table  VIII).  The  X^  and  X2  values  of  the  new  equation  differ 
very  little  from  the  first.  The  value  for  "a"  is  about  half,  since  100 
rather  than  200  8niq?le8  are  vred. )  When  this  equation  is  used  to  predict 
the  contour  counts  ^ich  ought  to  be  present,  and  the  prediction  compared, 
to  the  actual  number  of  contour  couTits,  a  further  coiparison  of  the  two 
halves  of  the  sanple  is  possible  (Table  IX). 

Table  IX  is  a  tabulation  .of  errors  by  masnitude,  whether  plus  or 
minixs.i  Because  of  the  rtlatively  wide  range  of  errors,  their  conparison 
is  best  made  by  the  use  of  a  cunnilatlve  curte  chart  of  errors.  Wiis  is 
done,  in  Figure  12,  based  on  the  data  from  Table  IX,  Bxet^t  at  the  ends 
of  the  distribution:,  the  distrlbutiosi  of  errors  is  similar  in,  both  grows- 
Underestimates  greater  than  25^  overestimates  greater  than  5r 

Involve  only  10,,  of  the  200  cases.  Either  the  symbol  indicating 


.“*cj  v:"' 


•rrom  in  Orpi^  cr  IImi  ispiilbeCL^'  *o**  ixM-ernii^  mrrorw  in  OroMfi 
couia  b«  u««d  ft«  fttidM  for  drtirtAf  •  mocih  Urn  te  rtf^rsMTit  ih»  4i*» 
triibiAtlan  of  orraro  In  •Itim’  or  boOi  mi^Ioo.  Ajk  iMu^d  bo  iai|>oeiod» 
tho  frootoot  mrrvra  oecttiTod  in  tho  ooao  coocr&phleoX  loeation  ai 
Um»  orrora  nopf  in  Fiforo  11.  Both  sM|>1o  Qrwfm  1  and  21  oahiblt 
tho  vum  diairiouiion  of  orrora  aial  aro  biaaod  towml  uodorpnaiictlon. 
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